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The incoherent regime of Multiple Andreev Reflections (MAR) is studied in long diffusive SNS 
junctions at applied voltages larger than the Thouless energy. Incoherent MAR is treated as a 
transport problem in energy space by means of a circuit theory for an equivalent electrical network. 
The current through NS interfaces is explained in terms of diffusion flows of electrons and holes 
through "tunnel" and "Andreev" resistors. These resistors in diffusive junctions play roles analogous 
to the normal and Andreev reflection coefficients in OTBK theory for ballistic junctions. The theory 
is applied to the subharmonic gap structure (SGS); simple analytical results are obtained for the 
distribution function and current spectral density for the limiting cases of resistive and transparent 
NS interfaces. In the general case, the exact solution is found in terms of chain-fractions, and 
the current is calculated numerically. SGS shows qualitatively different behavior for even and odd 
subharmonic numbers n — 2A/eV, and the maximum slopes of the differential resistance correspond 
to the gap subharmonics, eV = 2A/n. The influence of inelastic scattering on the subgap anomalies 
of the differential resistance is analyzed. 



I. INTRODUCTION 

The concept of multiple Andreev reflections (MARl 
was first introduced by Klapwijk, Blonder, and Tinkhamu 
in order to explain the subharmonic gap structure (SGS) 
on current-voltage characteristics of superconducting 
junctions. The theory was originally formulated for per- 
fect SNS junctions and then extended. to include the ef- 
fect of resistance of the SN interfaceo (OTBK theory). 
Within this approach, the subgap current transport is 
described in terms of ballistic propagation of quasiclas- 
sical electrons through the normal metal region, accom- 
panied by Andreev and normal reflections from specular 
NS boundaries. During every passage across the junc- 
tion, the electrons and the retro-reflected holes gain en- 
ergy equal to eV, which allows them eventually to escape 
from the SNS well. This energy gain results in strong 
quasiparticle nonequilibrium within the subgap energy 
region \E\ < A. 

OTBK theory gives a qualitatively adequate descrip- 
tion of dc current transport in voltage biased SNS junc- 
tions; however, its quantitative results have a rather lim- 
ited range of applicability. In short ballistic junctions 
with length d comparable with or smaller than the co- 
herence length (e.g., atomic-size junctions), the quantum 
coherence of subsequent Andreev reflections plays a cru- 



cial role leading to the ac Josephson effect. It has been 
shown that such a coherence also strongly modifies the 
dc current and SGSHu (coherent MAR regime). In fact, 
even in long ballistic SNS junctions (e.g. 2DEG-based 
devices), the coherence effects are important and give 
rise to resonant structures in the current due to Andreev 
quantization. In this respect, the quasiclassical OTBK 
theory, which does not include any coherence effects, may 
be qualified as a model for the incoherent MAR regime. 

One might expect that impurities could provide the 
conditions for incoherent MAR by washing out the An- 
dreev spectrum. However, this is not the case for a short 
diffusive junction, where appreciable Josephson coupling 
gives rise to coherent MAR.BQ The electron-hole coher- 
ence in the normal metal holds over a distance of the 
coherence length £e = y/tiD/2E from the superconduc- 
tor (D is the diffusion constant). The overlap of coherent 
proximity regions induced by both SN interfaces creates 
an energy gap in the electron spectrum of the normal 
metal, which plays the role of the level spacing in the 
ballistic case. In short junctions with a wide proximity 
gap of the order of the energy gap A in the supercon- 
ducting electrodes, the phase coherence covers the entire 
normal region. 

An incoherent MAR regime will occur in long diffusive 
SNS junctions with a small proximity gap of the order of 
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Thouless energy E Th - hV/d 2 < A.I3 If the applied volt- 
age is large on a scale of the Thouless energy, eV ^> i?Th, 
then the coherence length fg is much smaller than the 
junction length at all relevant energies E ~ min(eV, A). 
In this case, the proximity regions near the SN inter- 
faces become virtually decoupled and the Josephson os- 
cillations are strongly suppressed. At the same time, as 
soon as the inelastic mean free path exceeds the junction 
length, the subgap electrons must undergo many incoher- 
ent Andreev reflections before they enter the reservoir. 
We emphasize that such incoherency is provided by the 
small coherence length at large enough voltages, while 
the intrinsic dephasing length can be arbitrarily large. In 
order to describe such an incoherent MAR regime, one 
has to operate with the electron and hole diffusion flows 
across the junction rather than with ballistic quasiparti- 
cle trajectories, and to consider the Andreev reflections 
as the relationships between these diffusive flows. 

The first step in extending the OTBK approach to dif- 
fusive. SNS structures was taken by Volkov and Klap- 
wijk,Q who derived recurrence relations between the 
boundary values of the distribution functions. In that 
study, only a weak nonequihbrium was considered, which 
implies suppression of MAR by inelastic relaxation. In 
the present paper, we focus on the opposite case of strong 
nonequihbrium in the developed MAR regime, which re- 
sults in the appearance of SGS on I-V characteristics 
of the diffusive SNS junctions!! Following the interpreta- 
tion of MAR as a transport problem in energy spaceQ'E3, 
we analyze it by formulating an equivalent network in the 
spirit of Nazarov's circuit theoryO Within this approach, 
the energy-dependent tunnel and Andreev resistances of 
an equivalent circuit play roles similar to the normal and 
Andreev reflection probabilities in OTBK theory, and the 
effective voltage source is represented by Fermi reservoirs. 

The paper is organized as follows. In Section II, we de- 
rive the equations for incoherent MAR from the general 
Keldysh equations. In Sections III and IV, the circuit 
representation is formulated; some applications are con- 
sidered in Section V. The SGS in junctions with resistive 
interfaces is calculated in Section VI. The complete so- 
lution of the problem suitable for numerical calculation 
of the I-V characteristics is obtained in Section VII by 
using a chain-fraction techniques In Section VIII, we dis- 
cuss limitations on the MAR regime imposed by inelastic 
processes. 

II. MICROSCOPIC BACKGROUND 

The system under consideration consists of a normal 
channel (0 < x < d) confined between two voltage bi- 
ased superconducting electrodes, with the elastic mean 
free path I much shorter than any characteristic size 
of the problem. In this limit, the microscopic analy- 
sis of current transport can be performed within the 
framework of thej-diffusive equations of nonequihbrium 
superconductivityEJ for the 4x4 supermatrix Keldysh- 



Green function G(tit2,x): 

[H, G] = iKDdJ, J = Gd x G, G 2 = 1, (1) 

H = l[iha z d t - e<j>(t) + A(t)], A = Ae ia * x ia y , (2) 

where A is the modulus and x is the phase of the or- 
der parameter, and 4> is the electric potential. The Pauli 
matrices <7j operate in the Nambu space of 2 x 2 ma- 
trices denoted by "hats" , and the products of two-time 
functions are interpreted as their time convolutions. The 
junction length d is assumed to be smaller than the in- 
elastic and phase-breaking lengths, which allows us to 
exclude the inelastic collisions from our consideration at 
this stage; their role will be discussed later. The electric 
current / per unit area is expressed through the Keldysh 
component J K of the supermatrix current J: 

I(t) = (irhcr N /4:e)Tr<j z J K (tt,x), (3) 

where un is the conductivity of the normal metal. 

At the SN interface, the supermatrix G satisfies the 
boundary conditionllj 

(a N J) ±0 = (2R SN )- 1 [G^o,G +0 }, (4) 

where the indices ±0 denote the right and left sides of 
the interface and Rsn is the interface resistance per 
unit area in the normal state, which relates to, e.g., a 
Schottky barrier or mismatch between the Fermi veloci- 
ties. Within the model of infinitely narrow potential of 
the interface barrier, U(x) = H8(x), the interface resis- 
tance is related to the-barrier strength Z = H(Tivf)~ 1 
as Rsn — 21Z 2 /3(Tn£-3 It has been shown in Ref. [l5| 
that Eq. (|4|) is valid either for a completely transparent 
interface (Rsn — > 0, G+o = G_o) or for an opaque bar- 
rier whose resistance is much greater than the resistance 
R(l) — 1/un of a metal layer with the thickness formally 
equal to /. 

According to the definition of the supermatrix G, 

G= ( § %)iG K =g R f-fg A , (5) 

Eqs. (|]) and (|4|) represent a compact form of separate 
equations for the retarded and advanced Green's func- 
tions g R ' A and the distribution function f = f+ + cr z f-- 
Their time evolution is imposed by the Josephson rela- 
tion x{t) — 2eVt for the phase of the order parameter in 
the right electrode (we assume \ — in the left termi- 
nal). This implies that the function G(tit2,x) consists of 
a set of harmonics G(E n , E m , x), E n = E + neV, which 
interfere in time and produce the ac Josephson current. 
However, when the junction length d is much larger than 
the coherence length £g at all relevant energies E > eV, 
we may consider coherent quasiparticle states separately 
at both sides of the junction, neglecting their mutual in- 
terference and the ac Josephson effect. Thus, the Green's 
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function in the vicinity of left SN interface can be approx- 
imated by the solution g = a z cosh 9 + i<r y sinh 9 of the 
static Usadel equations for a semi-infinite SN structure,^ 
with the spectral angle 9(E, x) satisfying the equation 

ttmh[9{E,x)/A] = tanh[6» A r( J B)/4]cxp(-x/^ B Vi), (6) 
with the boundary condition 



where 



W^fiA/Esmh(9 N - 9 S ) + 2sinh(6V 2 ) = 0. (7) 

The indices S, N in these equations refer to the supercon- 
ducting and the normal side of the interface, respectively. 
The dimensionless parameter W in Eq. (0), 



W = 



Rsn 



dr 



Rsn 
Rn 



(8) 



where R^ = R(d) = d/apf is the resistance of the nor- 
mal channel per unit area, has the meaning of an effective 
barrier transmissivity for the spectral functions .113 Note 
that even at large barrier strength Z » 1 ensuring the 
validity of the boundary conditions Eq. @)fB the effec- 
tive transmissivity W ~ (£a/1)Z~ 2 of the barrier in a 
"dirty" system, I -C £a, could be large. In this case, the 
spectral functions are virtually insensitive to the pres- 
ence of a barrier and, therefore, the boundary conditions 
Eqs. (0) can be applied to an arbitrary interface if we 
approximately consider high-transmissive interfaces with 
W > £aA ^ 1 as completely transparent, W = oo. 
For low transmissivity, IF <§; 1, Eq. (^) can be analyzed 
within a perturbative approach (see the Appendix). At 
arbitrary W, Eq. (Q) should be solved numerically. 

The distribution functions f±(E,x) are to be consid- 
ered as global quantities within the whole normal channel 
determined by the diffusive kinetic equations 



d x [D±(E,x)d x f±(E,x)] =0, 
with dimensionless diffusion coefficients 

D+ = (1/4) Tr(l - g R g A ) = cos 2 ho. 9, 



D- = (l/4)Tr(l 



a z g R a z g A ) 



cosh Re 9. 



(9) 

(10a) 
(10b) 



Assuming the normal conductance of electrodes to be 
much greater than the junction conductance, we con- 
sider them as equilibrium reservoirs with unperturbed 
spectral characteristics, 9s = Arctanh(A/i?), and equi- 
librium quasiparticle distribution, fs{E) = fo{E) = 
tanh(_B/2T). Within this approximation, the boundary 
conditions for the distribution functions in Eq. (|^) at 
x = read 

a N D + d x f+(E,0) = G + (E)[f + (E,0)-f (E)], (11) 



* N D-d x f-(E,0) = G-(E)f-(E,0), 



(12) 



G ± (E) 



p— i 



(N s N n tM±M±) 



(13) 



N(E) = Re(cosh0), M + {E) + iM~{E) = sinh.9. (14) 

At large energies, \E\ 3> A, when the normalized den- 
sity of states N(E) approaches unity and the condensate 
spectral functions M ± (E) turn to zero at both sides of 
the interface, the conductances G±(E) coincide with the 
normal barrier conductance; within the subgap region 
\E\ < A, G+(E) = 0. 

Similar considerations are valid for the right NS in- 
terface, if we eliminate the explicit time dependence of 
the order parameter in Eq. (Q), along with the potential 
of right superconducting electrode, by means of a gauge 
transformationt3 

G(t\t2, x) = exp(ia z eVti)G(tit2, x) exp(— icr z eVt2). (15) 

As a result, we arrive at the same static equations and 
boundary conditions, Eqs. (p|)-(fl4|), with x — * d — x, 

for the gauge-transformed functions g(E, x) and /(£", x). 
Thus, to obtain a complete solution for the distribution 
function /_, which determines the dissipative current 



2e 



dED_d x f_, 



(16) 



we must solve the boundary problem for f(E,x) at the 
left SN interface, and a similar boundary problem for 

f(E, x) at the right interface, and then match the dis- 
tribution function asymptotics deep inside the normal 
region by making use of the relationship following from 
Eqs. ©, ©: 



f(E,x)=f(E + a z eV,x) 



(17) 



III. 



CIRCUIT REPRESENTATION OF 
BOUNDARY CONDITIONS 



In order for this kinetic scheme to conform to the con- 
ventional physical interpretation of Andreev reflection in 
terms of electrons and holes, we introduce the following 
parameterization of the matrix distribution function, 



f(E,x) = l- 



n e (E,x) 






i h (E,x) 



(18) 



where n e and n h will be considered as the electron and 
hole population numbers. Deep inside the normal metal 
region, they acquire rigorous meaning of distribution 
functions of electrons and holes. In equilibrium, the func- 
tions n e ' h approach the Fermi distribution. In this rep- 
resentation, Eqs. (H) take the form 
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D±{E,x)d x n±{E,x) = const = -I ± (E)/a N , (19) 

where n± = n e ± n , and they may be interpreted as 
conservation equations for the (specifically normalized) 
net probability current I + of electrons and holes, and 
for the electron-hole imbalance current i_. Furthermore, 
the probability currents of electrons and holes, defined 
as I e ' h — (l/2)(/ + ± 7_), separately obey the conserva- 
tion equations. The probability currents I e,h are nat- 
urally related to the electron and hole diffusion flows, 
I e ' h = —aNd x n e ' h , at large distances x ^> £e from the 
SN boundary. Within the proximity region, x < £e, each 
current I e h generally consists of a combination of both 
the electron and hole diffusion flows, 

I e > h = -(<T N /2)[(D + ±D-)d x n e + (D + TD-)d x n h ] , (20) 

which reflects coherent mixing of normal electron and 
hole states in this region. 




b) c) 

FIG. 1. Elementary equivalent circuits representing boun- 
dary conditions Eq. ( J2l| ) for the electron and hole population 
numbers n e ' h (E,0) and probability currents I e ' h (E), at en- 
ergies outside the gap, \E\ > A (a), and within the subgap 
region, \E\ < A (b); equivalent network in (x,i?)-space for in- 
coherent MAR in SNS junction (c). Filled and empty symbols 
stand for electron- and hole-related elements, respectively; 
half-filled squares denote Andreev resistors; U„ = np{E n ). 

In terms of electrons and holes, the boundary condi- 
tions in Eqs. (|ll|), @ read 

I e ' h = G T {n F - rf- h ) T G A {n e - n h ), (21) 

where 

G T = G+, G a = (G_ - G+)/2. (22) 



Each of the equations Eq. ( pl| ) may be clearly in- 
terpreted as a Kirchhoff rule for the electron or hole 
probability current flowing through the effective circuit 
(tripole) shown in Fig. 0(a). Within such an interpre- 
tation, the nonequilibrium populations of electrons and 
holes n e,h at the interface correspond to "potentials" of 
nodes attached to the "voltage source" - the Fermi dis- 
tribution np(E) in the superconducting reservoir - by 
"tunnel resistors" R T (E) = G^ T 1 (E). The "Andreev re- 
sistor" Ra(E) = G^ 1 (E) between the nodes provides 
electron-hale conversion (Andreev reflection) at the SN 
interfaceO 

The circuit representation of the diffusive SN inter- 
face is analogous to the scattering description of ballis- 
tic SN interfaces: the tunnel and Andreev resistances in 
the diffusive case play the same role as the normal, and 
Andreev reflection coefficients in the ballistic case.EJ For 
instance, for \E\ > A [Fig. |l](a)], the probability current 
I e is contributed by equilibrium electrons incoming from 
the superconductor through the tunnel resistor and 
also by the current flowing through the Andreev resistor 
Ra as the result of hole-electron conversion. Within the 
subgap region, \E\ < A, [Fig. 0(b)], the quasiparticles 
cannot penetrate into the superconductor, Rt = oo, and 
the voltage source is disconnected, which results in de- 
tailed balance between the electron and hole probability 
currents, I e = — I h (complete reflection). For the per- 
fect interface, Ra turns to zero, and the electron and 
hole population numbers become equal, n e = n h (com- 
plete Andreev reflection) . The nonzero value of the An- 
dreev resistance for R$n ^ accounts for suppression 
of Andreev reflection due to the normal reflection by the 
interface. 

Detailed information about the boundary resistances 
can be obtained from asymptotic expressions for the 
bare interface resistances R±(E) = G^. 1 (E) (see the Ap- 
pendix) and numerical plots of R±(E) in Fig. || In par- 
ticular, R±(E) turns to zero at the gap edges due to 
the singularity in the density of states which enhances 
the tunneling probability. Furthermore, the imbalance 
resistance R-{E) approaches the normal value Rsn at 
E — > due to the enhancement of the Andreev reflection 
at small energies, which results from multiple coherent 
backscattering of quasiparticles by the impurities within 
the proximity region. This property is the reason for the 
re-entrant behavior of the conductance of high-resistive 
SIN systemstlcj at low voltages. In the MAR regime, one 
cannot expect any reentrance since quasiparticles at all 
subgap energies participate in the charge transport even 
at small applied voltage. 

The diffusion coefficients D± in Eq. ([[o]) turn to unity 
far from the SN boundary, and therefore the population 
numbers n e ' h become linear functions of x, 

n e ' h {E, x) w n e ' h (E, 0) - R N I e ' h (E)x/d. (23) 

This equation defines the renormalized population 
numbers n e,h (E, 0) at the NS interface, which differ from 
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n e ' (E, 0) due to the proximity effect, as shown in Fig. 
||. These quantities have the meaning of the true elec- 
tron/hole populations which would appear at the NS in- 
terface if the proximity effect had been switched off. ft 
ispossible to formulate the boundary conditions in Eq. 
( |2f| ) in terms of these population numbers by including 





FIG. 2. Energy dependencies of bare interface resistances 
R±(E) = Gjjj (E) (a,b), bare Andreev resistance Ra(E) (c) 
and normalized proximity corrections fi± (E) [insets in (a) and 
(b)], for different values of the resistance ratio r = Rsn/Rn 
and d/^A = 5. 



the proximity effect into renormalization of the tunnel 
and Andreev resistances. To this end, we will associate 
the node potentials with renormalized boundary values 
n e - h (E,0) = (l/2)[n+(E,Q) ±n-(E,0)] of the popula- 
tion numbers, where n±(E, 0) are found from the exact 
solutions of Eqs. (|l9|), 



n±(E,0) = n ± (E,0) - m±{E)I ± (E). 



(24) 



Here m± (E) are the proximity corrections to the normal 
metal resistance at given energy for the probability and 
imbalance currents, respectively, 



m±(E) = ±R N (U/d)n±(E), 



M±(£)= / ^\Dl\E,x) 

JO ?A 



1 > 0, 



(25a) 



(25b) 



see the insets in Fig. |](a,b). It follows from Eq. (|2~i| ) 
that the same Kirchhoff rules as in Eqs. (|2|), ( p2| ) hold 

for n e,h (E,0) and I e - h (E), if the bare resistances R± are 
substituted by the renormalized ones, 

R±(E) — > R~±(E) = R±(E) + m±(E). (26) 




n(E,0) 
n(E,0) 



FIG. 3. Qualitative behavior of population numbers within 
the normal channel (solid curve). The edge distortions of 
the linear ^-dependence of population numbers, Eq. (^), oc- 
cur within the proximity regions. The difference between the 
boundary population numbers n(E,0), n(E,d) and their ef- 
fective values n(E,0), n(E,d) for true normal electrons and 
holes is included in the renormalization of the boundary re- 
sistances, Eq. (ptj). 

The energy dependence of the renormalized boundary 
resistances Rt(E) and Ra(E) is illustrated in Fig. || In 
some cases, there is an essential difference between the 
bare and renormalized resistances, which leads to qualita- 
tively different properties of the SN interface for normal 
electrons and holes compared to the properties of the bare 
boundary. Let us first discuss a perfect SN interface with 
Rsn 0. Within the subgap region \E\ < A, the bare 
tunnel resistance Rt is infinite whereas the bare Andreev 
resistance Ra turns to zero; this corresponds to com- 
plete Andreev reflection, as already explained. However, 
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the Andreev resistance for normal electrons and holes, 
Ra{E) — 2m^(E), is finite and negativepj which leads 
to enhancement of thejprmal metal conductivity within 
the proximity region.EH At \E\ > A, the bare tunnel 
resistance Rt is zero, while the renormalized tunnel resis- 
tance Rt(E) = m + (E) is finite (though rapidly decreas- 
ing at large energies). This leads to suppression of the 
probability currents of normal electrons and holes within 
the proximity region, which is to be attributed to the ap- 
pearance of Andreev reflection. Such a suppression is a 
global property of the proximity region in the presence of 
sharp spatial variation of the order parameter, and it is 
similar to the over-the-barrier Andreev reflection in the 
ballistic systems. In the presence of normal scattering 
at the SN interface, the overall picture depends on the 
interplay between the bare interface resistances R± and 
the proximity corrections m± ; for example, the renormal- 
ized tunnel resistance Rt(E) diverges at \E\ — ► A, along 
with the proximity correction m + (E), in contrast to the 
bare tunnel resistance Rt(E). This indicates complete 
Andreev reflection at the gap edge independently of the 
transparency of the barrier, which is similar to the sit- 
uation in the ballistic systems where the probability of 
Andreev reflection at \E\ = A is always equal to unity. 




lof 



r = 1 
' 0.2 


' ' / 




i.i, 



b) 



E/A 



FIG. 4. Energy dependencies of the renormalized interface 
resistances: tunnel (probability) resistance Rt(E) (a) and 
Andreev resistance Ra(E) (b), for d = 5£a- 



IV. ASSEMBLING MAR NETWORKS 

To complete the definition of an equivalent MAR net- 
work, we have to construct a similar tripole for the right 
NS interface and to connect boundary values of popula- 
tion numbers (node potentials) using the matching con- 
dition Eq. ( [T~T| ) expressed in terms of electrons and holes: 



l (E,x) = ri e ' h (E±eV,x). 



(27) 



Since the gauge-transformed distribution functions f± 
obey the same equations Eq. (p|)-(p^), the results of the 
previous Section can be applied to the functions n e ' h (E) 
and —I e ' h (E) (the minus sign implies that / is associated 
with the current incoming to the right-boundary tripole). 
In particular, the asymptotics of the gauge-transformed 
population numbers far from the right interface are given 
by the equation 



(E,d)+R N I e < h (E)(l-x/d). (28) 



l (E,x) 



After matching the asymptotics in Eqs. (|23j) and (|2 
by means of Eq. (pTl), we find the following relations 



P' h (E) = P' h (E±eV), 



n e ' h 



(E,0) 



-e.h 



{E±eV,d) =R N I e > h (E) 



(29) 



(30) 



From the viewpoint of the circuit theory, Eq. ( ppj ) may 
be interpreted as Ohm's law for the resistors Rn which 
connect energy-shifted boundary tripoles, separately for 
the electrons and holes, as shown in Fig. |l|(c). 

The final step which essentially simplifies the analysis 
of the MAR network, is based on the following observa- 
tion. The spectral probability currents I e,h yield opposite 
contributions to the electric current in Eq. (|l6|), 



1 r°° 



(31) 



due to the opposite charge of electrons and holes. At the 
same time, these currents, referred to the energy axis, 
transfer the charge in the same direction, viz., from bot- 
tom to top of Fig. |l](c), according to our choice of pos- 
itive eV. Thus, by introducing the notation I n (E) for 
an electric current entering the node n with the energy 
E n = E + neV, as shown by arrows in Fig. 0(c), 



In{E) 



/ e (^„-l), 
-l' l (En), 



--2k + l, 

= 2k, 



(32) 



we arrive at an "electrical engineering" problem of cur- 
rent distribution in an equivalent network in energy space 
plotted in Fig. ||, where the difference between electrons 
and holes becomes unessential. The bold curve in Fig. |B| 
represents a distributed voltage source - the Fermi dis- 
tribution uf(E) connected periodically with the network 
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nodes. Within the gap, |2? n | < A, the nodes are discon- 
nected from the Fermi reservoir and therefore all partial 
currents associated with the subgap nodes are equal. 

Since all resistances and potentials of this network de- 
pend on E n — E + neV, the partial currents obey the 
relationship I n (E) = Ik[E + (n — fc)eV] which allows us 
to express the physical electric current, Eq. (|3l]), through 
the sum of all partial currents /„ flowing through the nor- 
mal resistors Rn, integrated over an elementary energy 
interval < E < eV: 

I=— dE[h{E) + I (E)} = - / dEJ(E), (33) 



J(E)= J2 UE)- 



(34) 



The spectral density J(E) is periodic in E with the pe- 
riod eV and symmetric, J(—E) = J(E), which follows 
from the symmetry of all resistances with respect to E. 




FIG. 5. MAR network of Fig. [j](c) in energy space. The 
nodes outside the gap are connected with the distributed volt- 
age source n F (E) (bold curve); the subgap nodes are discon- 
nected from the voltage source. 



As soon as the partial currents are found, the popu- 
lation numbers can be recovered by virtue of Eqs. (19) 
©, ©, and ©: 

n e ' h (E,x) 



n e ' h 



(E,0)TR N h,o{E)x/d, 



(35) 



n 



1 



(E, 0) = n F - - [R+ (h -h)±R- (Ji + Jo)] (36) 



at \E\ > A. Within the subgap region, Eq. (J36p is inap- 
plicable due to the undeterminacy of product R + (1% —Jo). 
In this case, one may consider the subgap part of the net- 
work as a voltage divider between the nodes nearest to 
the gap edges, having the numbers — iV_, N + , respec- 
tively, where N± = Int[(A =p E)/eV] + 1, Int(x) denoting 
the integer part of x. Then the boundary populations at 
\E\ < A become 



r?> h (E,0) = n L ' R (E ±N± ) 
± Jo 



N±-l 

N ± R N + RA(E± k ) 
fe=i 



(37) 



where R, L indicate the right (left) node of the tripole, 
irrespectively of whether it relates to the left (even n) 
or right (odd n) interface. The physical meaning of 



,R,L 



(E n ), however, depends on the parity of n: 



(En) 



I n e > h (E n ,0), n = 2k, 
' [ n H,e (E n ,d), n = 2k + l. 

in Eq. ( |37| ) can be found from Eq. ( |3^ ) 
which is generalized for any tripole of the network in Fig. 
outside the gap as 



The values 



,R,L 



(E n ) = n F (E n ) 



(1/2) [R + (E n )(I n+1 -I n ) 
±R^{E n ){I n+l +/„)] 



(39) 



The circuit formalism can be simply generalized to the 
case of different transparencies of NS interfaces, as well as 
to different values of A in the electrodes. In this case, the 
network resistances become dependent not only on E n 
but also on the parity of n. As a result, the periodicity of 
the current spectral density doubles: J(E) = J(E+2eV), 
and, therefore, J(E) is to be integrated in Eq. ( |33| ) over 
the period 2eV, with an additional factor 1/2. 



V. SIMPLE APPLICATIONS 

A helpful example of an asymmetric junction which al- 
lows us immediately to obtain an analytical solution is 
given by the SNN structure. The problem of calculation 
of its conductivity is inherently static and therefore may 
be completely solved for any junction length. If the lat- 
ter is much larger than the coherence length, the circuit 
approach of the previous Section is applicable without 
restrictions. Due to the absence of superconducting cor- 
relations at the right NN interface, odd Andreev resistors 
are eliminated and, therefore, the whole network may be 
split into separate finite circuits located around even (su- 
perconducting) nodes, as shown in Fig. ^(a); moreover, 
odd tunnel resistances are to be considered as normal 
ones. After some simple algebra, we obtain the sum of 
partial currents in each subcircuit, 



12k 



l2k+l — 



n F (E 2 k- 



nF(E 2 k+i) 



Rn + Rsn + R-{E: 



(40) 



2 A- ) 



which leads to a well known formula for the I-V charac- 
teristics of a long SNN junction:!! 



1 

~2e 



2eV 



dE 



E 

k—~ oo 



n F {E 2 k-i) ~ n F {E; 



2fc+l J 



dE 



Rn + Rsn + R-(E 2k ) 
n F (E - eV) -n F (E + eV) 
R N + Rsn + R-{E) 



(41) 



If the junction is short enough (d and £e are compa- 
rable), one might naively expect some kind of proximity- 
induced Andreev scattering at the right NN interface, fol- 
lowed by MAR and SGS anomalies in the I-V character- 
istic. However, the circuit theory rejects this assumption 
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at once: since the condensate spectral functions M ± , Eq. 
(|l4|), disappear in the normal terminal, the conductivi- 
ties G± become equal, and the Andreev channel becomes 
closed (Ga — 0) at the NN interface at any length of the 
junction. Thus, the circuit model of charge transport in 
Fig. |6|(a) remains valid, with a few modifications: (i) the 
spectral angle 8(E, x) has to be found from the Usadel 
equation for the finite interval < x < d, (ii) the prox- 
imity corrections m± at the left SN interface are to be 
expressed through the integrals 



m±(E) 



±R N J ^\Dl\E : x)-l\, 



(42) 



instead of Eq. (|25|), and (iii) odd tunnel resistances Rsn 
should be replaced by the energy-dependent bare resis- 
tances R+(E) [or, equivalently, R_(E)]. From this point 
of view, the entire channel represents a global "Andreev 
reflector" for normal electrons and holes incoming from 
the right reservoir. 



[l + njr(£i)-n F (£-i)]ir 



(43) 



where R = Rn + 2Rsn is the net normal resistance of 
the junction. In summary, we obtain another well-known 
result E9 



^ dE n F (E-eV)-n P (E + eV) 
Rn + Rsn + R-(E) 



eR Jo 

-R_{E)]* A f°° dE Xm-X-W (44) 
eRJo R N +R SN + R-(E) 

It is of interest to note that the net resistance Rt = 
i?+ for the probability current never enters final results 
in these examples and, therefore, the superconducting 
modification involves only the imbalance resistance R_. 
In other words, only the evolution of the imbalance n_ 
between the electron and hole populations is relevant for 
the charge transport in such cases. 



+ [ hk Y ! 2k+<l \ + 



VI. MAR TRANSPORT 



U 2k-1 U 2k U 2k+1 
a) 

b) 




+ . 
U 2 




FIG. 6. Simplified circuits used for calculation of (a) I-V 
characteristic of SNN junction; (b) excess (deficit) current 
in SNS junction; (c) I-V characteristic of SNS junction with 
low-transparent interfaces. 

The next simple application of this circuitry is given 
by calculation of the excess (deficit) current in SNS junc- 
tions, i.e., the difference I ex = I(V) — V/R between the 
currents in superconducting and normal state at large 
voltages V ^> A/e. Assuming the integration in Eq. j3^ ) 
to be reduced to the interval < E < eV/2 by making 
use of the symmetry J(E) — J(—E) — J(eV — E), we 
note that the Andreev conductances are negligibly small 
for all nodes with n ^ (E n ^> A). Thus, the circuit in 
Fig. [| may be split, as in the case of the SNN junction 
discussed above, into the chain of separate circuits shown 
in Fig. ^|(b). The contribution of the central circuit is de- 
scribed by Eq. (40) with k = 0, whereas the other parts 
are to be considered as normal circuits and represent con- 
tribution of thermally excited quasiparticles: 



Proceeding with the analysis of current transport 
through the SNS junction at arbitrary voltages, we first 
discuss the case of low-transparent barriers, W <C 1. 
We note that in practice this case is relevant for a wide 
range of junctions both with high interface resistance, 
Rsn S> Rn, and comparatively low interface resistance 
Rsn "C Rn- Indeed, according to Eq. (|J), the ra- 
tio Rn/Rsn = Wd/^A, being proportional to W, con- 
tains also the large parameter d/^A- Therefore, the limit 
W "€L 1 covers most of the practically interesting situa- 
tions, < Rn/Rsn *C d/^A, and only the case of very 
small interface resistances, Rn/Rsn > d/^A ^> 1, re- 
quires special consideration. 

At W <C 1, the proximity effect is essentially sup- 
pressed and the calculations can be performed on the ba- 
sis of a simplified model of the equivalent network, which 
nevertheless provides a quantitative description of I(V). 
Due to the sharp increase in the Andreev resistance at 
\E\ > A [see Fig. §(b)], all Andreev resistors outside the 
gap can be excluded, and we arrive at the sequence of 
subcircuits shown in Fig. ^|(c). The central circuit be- 
tween the nodes — 7V_ and N + includes N + + N- normal 
and N + + iV_ — 1 Andreev resistors within the gap, as 
well as two tunnel resistors at the circuit edges. The total 
resistance Ra of this circuit is 

Ra(E) = (JV+ + N_)R N + E (E k ) 

-JV_<fe<7V + 

+Rt(E n+ ) +Rt(E_ n _), (45) 

and the current Iq through this circuit is given by Ohm's 
law: 



I {E) = [n F {E- N -) - n F (E N+ )]/RA(E). 



(46) 



8 



Thus, the contribution of this circuit to the current spec- 
tral density, Eq. (fj), is (JV+ +iV_)J . 

The current of thermal excitations is carried by the 
side circuits (n > N + , n < — iV_): 

n F (E n -i) - n F (E n ) 



(47) 



R N + R T (E n ) + R T (E n _x) 
From Eq. (M9) we obtain a simple formula at T < A: 



I = 



dE 



, Rmar(E) 



Ra(E) 
N + + N_' 



(48) 



where Rmar(E) has the meaning of the effective resis- 
tance of the subgap region for the physical electric cur- 
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q 2MeV 

FIG. 7. I-V characteristics (a) and differential resistance 
vs inverse voltage (b) of SNS junctions with low-transparent 
interfaces, W — 0.1 and W = 0.2, at g(/£a = 5. Dashed lines 
denote low- voltage asymptotics of the I-V curves, Eq. (|49j). 



parameter W was chosen to be equal to 0.1 and 0.2 
at c?/£a — 5, which corresponds to the resistance ratio 
r = Rsn/ Rn equal to 2 and 1, respectively. In our calcu- 
lation of Rt : a{E) in Eq. ([45|), we used the asymptotic Eq. 
( |A3| ) for the bare resistances Rt,a (E) at W <C 1 , neglect- 
ing small proximity corrections m±(E) ~ Rn(£,a/ 'd)W 2 , 
Eq. ( |A2bp . The results are in good agreement with those 
obtained on the basis of exact calculations [see further 
Eq. (p3|)]. The smeared steps in the I-V characteristic in- 
dicate steplike increase in the number of subgap Andreev 
resistors in Eq. ([45|). The sharp peaks (dips) in dV/dl 
arise from the rapidly varying contribution of the nodes 
which simultaneously cross the gap edges, and therefore 
both the edge Andreev resistances undergo strong sup- 
pression. A certain contribution also comes from the edge 
tunnel resistances which also show singular behavior at 
— > A. The peaks are more pronounced for even sub- 
harmonics, when the middle Andreev resistor crosses the 
Fermi level, E = 0, and its value is suppressed simultane- 
ously with the edge Andreev resistors. Careful analysis 
shows that at the gap subharmonics, eV n — 2A/n, the 
second derivative d 2 V/dI 2 has sharp maxima. 

The magnitude of dV/dl strongly depends on the num- 
ber of large Andreev resistors which contribute to -Rmar- 
At eV < A, at least one Andreev resistor appears far 
from the "resonant" points E — 0, ±A where Ra sharply 
decreases. Thus, the net subgap resistance Rmar{E) 
remains large (~ Rsn /W) at any energy, which re- 
sults in large differential resistance dV/dl ~ Rsn/W 
at these voltages. In the vicinity of the second sub- 
harmonic, eV ~ A, the current transport involves a 
high-transmissive circuit with three Andreev resistors 
located near the resonant points, which yields a much 
smaller value of dV/dl ~ Rsn- The same effect oc- 
curs at eV ~ 2A when the resonant circuit contains 
two suppressed Andreev resistances at the gap edges. At 
eV > 2A the differential resistance is basically deter- 
mined by quasiparticles which overcome the energy gap 
without Andreev reflections, and it turns to the normal 
value R at large voltages. 

At low voltage, the amplitude of the oscillations of the 
differential resistance decreases and the asymptotic value 
of dV/dl at V <§; A/e can be found analytically from 
Eqs. ([45|), (f48|), by replacing the sum in Eq. ( |45| ) with an 
energy-independent integral. As a result we get 



dV/dI{Q)^R N - 



1 f A 

— J dER A {E)=R N - 



16a/2 
'21W 



Rsn- (49) 



Since each circuit in Fig. ^(c) represents a separate 
voltage divider, we easily obtain the boundary values of 
the population numbers. If the node n = belongs to 
the central circuit (—A — eV < E < A for electrons and 
— A < E < A + eV for holes), we have 



In Fig. ^ we present the I-V characteristics and the 
differential resistance vs inverse voltage at zero temper- 
ature, calculated numerically by means of Eq. ([48|). The 



—e.h 



(E,0) = n F (E± N± )± 



n F {E. 



n - 



n F (E N+ ) 



Ra(E) 
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Rt(E± n± ) + N±R 



N 



N±-l 

fe=l 



(50) 



and, in the opposite case, 

R T (E)[n F (E)~n F (E±eV)] 



n e - h 



(E, 0)-n F (E)- Rn + RHE) + Rt{e ± eV) 



As follows from Eqs. (p0|), (|51|), at low temperatures, 
the energy distribution of quasiparticles within the re- 
gion —A — eV < E < A + eV has a steplike form (Fig. 
|8|) , which is qualitatively similar to, but quantatively dif- 
ferent from, that found in OTBK theory.B The number 
of steps increases at low voltage, and the shape of the 
distribution function becomes resemblant to a "hot elec- 
tron" distribution with the effective temperature of the 
order of A. This distribution is modulated due to the 
discrete nature of the heating mechanism of MAR, which 
transfers the energy from an external voltage source to 
the quasiparticles by energy quanta eV. Since the subgap 
probability current Iq, Eq. (ff6|), is strongly suppressed by 
large subgap resistance Ra ~ Rsn(N + + NJ)/W 3> Rn, 
the spatial variations of the population numbers, Eqs. 
(||), (§|), are negligibly small: n e ' h (E, x) - n e ' h (E, 0) ~ 
Rn/Ra < 1- 
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FIG. 8. Energy dependence of the electron population 
number n e (E, 0) at the left interface of the SNS junction with 
Rsn = Rn and d = 5£a, at voltage V = A/e and T = 0. 



VII. EXACT SOLUTION 

In the case of high transmittance of the NS interface, 
the partial currents outside the gap noticeably contribute 
to the net electric current even at T = 0. In this case, 
the I-V curves should be calculated on the basis of the 
exact solution of the recurrence relations between partial 
currents, 

In+irn + In-ir n -i-In(Pn + rn-i+r„) = U n -U n -i, (52) 

following from the Kirchhoff rules for an infinite network 
in Fig. @. Here U n (E) = n F {E n ) and 



Pn (E) = R N + R_ (£„_! )+R_(E n ), (53a) 

r n (E) = \R+(E n ) - R_{E n )}/2. (53b) 

By analogy with differential equations, we introduce 
the following ansatz: 

I n {E) = A+(E)I+(E) + A-{E)I-{E), (54) 

where 1^ are the fundamental solutions of the corre- 
sponding uniform equation, decreasing at n — ► ±oo, re- 
spectively, 



Sn-i,Q, n > 0, 
Szl , n<0, 



(55a) 



™={?- 1 :;°;^o; ^ 



Iq = I Q =1. The quantitic 



S mn (E) = H Sj (E), P mn (E) = J] Pj( E ) ( 56 ) 



are expressed through the products of chain fractions 

i+l/^n i Pn = I n / I n + 



s n = l£ + i/Ini Vn = I n l^n+ii defined by the recurrences 



Pn = 



(57a) 



Pn+i + r n + a n+ i p n + r„ + 6„-i 

a n = r n (l - s n ), b n = r„(l - p n ), (57b) 

with the boundary conditions s n ^ +OQ — > 0, p„^-oo —> 0. 
Within the gap, \E n \ < A, where r n — > oo, the values s n , 
p n are equal to 1, in accordance with the conservation of 
the subgap currents mentioned above. 

The coefficients A^ in Eq. (p4f) satisfy an equation fol- 
lowing from Eqs. (|2(), (|I|), (jf|), 

{(P ^jt+i - r„s-MA+]5 n , + [(pn+i 

-r n p n )SAZ +1 - r n p n 6A~]P~o = U n+1 - U n , (58) 

for n > (and similar for n < 0), where SA n — A n+ i — 
A n . The requirement of cancellation of terms with 5A n z +1 
in Eq. ( |58| ) allows us to completely determine A^. This 
yields first-order recurrences for A^ which lead to the 
formula for partial currents at n > 0, 

n 

In = A+I+ + A -/- + £ j k (S n -i,k ~ P-\ k ) , (59) 
fe=l 

j n {E) = ([/„_! - U n )/( Pn + a n + &„_!). (60) 

The undeterminacy of a n , in Eq. ( |60| ) within the 
subgap region [where s n = p n = 1 and r n — > oo, see Eq. 
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( |57b| )], is resolved by the recurrences a n = p n +i + Q-n+i, 
b n = Pn + &n-i following from Eqs. (jo"^ ) for \E n \ < A. 
These recurrences are to be continued until the nodes 
N + and — A_, respectively, are reached. As a result, we 
obtain a convenient representation for j n at (— N- + 1) < 
n < N+: 



UE) 



U n -\ — U n 



Ra ~ r N+ (l + s N+ ) - r_Ar_(l +P-N-) ' 



(61) 



The effective subgap resistance determined by the de- 
nominator in Eq. (JsTJ) differs from Ra in Eq. ( fl5| ) by extra 
terms describing leakage of the subgap current through 
the Andreev resistors outside the gap. 

The coefficients 
in Eq. (|59| ) at all n; for instance, for n 



Aq have to provide finite values of /„ 
-oo, the diver- 



gent products of p- in I and P 



each other: A Q = Y. k =o3kP} 



\ k should compensate 

'k-i,o- A similar procedure 
for negative n determines the value of Aq and results in 
the final formula for the partial current with arbitrary 
index n, 

I n (E)=j n (E)+ 3k{E)P k -l,n(E) 

k=n+l 
n-1 

+ J2 jk(E)S n ^, k (E). (62) 

k— — oo 

By making use of the relation p_ n (E) — s n (—E) fol- 
lowing from Eq. ( |57|) and taking into account the sym- 
metry of all resistances with respect to E, we obtain the 
net current spectral density, 



+ 30 



J ( E )= E \]n(E) + J2lME)S k , n (E) 



+j- n (E)S k+1 . n+1 {-E)} 



(63) 



equal to the junction length; the variation of is not 
critical for the fine structure of dV/dl. 




b) 
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FIG. 9. I-V characteristics (a) and differential resistances 
vs inverse voltage (b) of an SNS junction with high-transpa- 
rent interfaces. 



At low temperatures T <C A, only the term with n 
contributes to the sum in Eq. (|63|), 



J(E) = j (E) j 1 + E i S k-ho(E) + S k ,i (-E)] I . (64) 

Figure ^| shows the results of a numerical calcula- 
tion of I(V) and dV/dl for an SNS junction with high- 
transmissive interfaces, W = 1 and W — > oo, at zero tem- 
perature. In practice, due to rapid decrease of the coeffi- 
cients r n (E) in Eq. ([33]) at large energies, it is enough to 
calculate recurrences in Eq. ( p7|) starting from the max- 
imum number n max = iV_ + 2 and assuming the chain 

,.v. To 



fractions to be truncated, s r , 



0, at n > 



avoid formal singularity in m_(_B) at E — > 0, we intro- 
duce a small dephasing factor T which provides a cutoff 
of the coherence length £e — > y/KD/2{E + iT). The cor- 
responding dephasing length 1$ = y / HV/2T was chosen 



Similar to the case of low barrier transmittance, the 
current transport through an SNS junction with nearly 
perfect interfaces can be qualitatively explained within 
a simplified model of MAR, where the over-the-barrier 
(\E\ > A) Andreev reflection is ignored. Indeed, as fol- 
lows from Fig. ||, at Rsn the tunnel resistances 
Rt(E) outside the gap are much smaller than the An- 
dreev and/or normal resistances, except for a narrow re- 
gion around the gap edges, where Rt{E) diverges due to 
complete Andreev reflection. This allows us to assume all 
the normal resistors at \E\ > A to be connected directly 
to the "voltage source" np(E) and therefore to exclude 
all Andreev resistors in Fig. || outside the gap. As a 
result, we arrive at the sequence of subcircuits shown in 
Fig. ||(c), with R T = for side circuits. Thus, at T < A, 
the subgap current may be approximately described by 
Eqs. (f45|), (48), with the tunnel and Andreev resistances 
renormalized by the proximity effect. 
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Within this model, the SGS oscillations in the differ- 
ential resistance in Fig. ^ can be explained in the follow- 
ing way. As the voltage decreases, the subgap current, 
which approximately follows Ohm's law, undergoes an 
additional suppression in the vicinity of the gap subhar- 
monics, due to the presence of a high-resistive circuit 
with two large tunnel resistors located just at the gap 
edges. These current steps, being almost invisible in the 
I-V characteristic, manifest themselves as sharp dips in 
dV/dl. At even subharmonics, this effect is partially 
compensated by the middle negative Andreev resistor, 
which rapidly reduces the effective normal metal resis- 
tance due to the increase in the size of the proximity 
region at small energies. As a result, the even dips be- 
come less pronounced and, as long as the interface trans- 
parency increases, turn into small peaks. At low voltages, 
the differential resistance approaches a constant value, 
which can be estimated for perfect NS interface by the 
following expression similar to Eq. (|^|): 



dV/dI(0) « R N 



-A 



dE 
2A 



R N (1 - 2.64&A0 , R S n = 0. 



(65) 



Unlike the ballistic SNS n junction,El'H but similar to 
short diffusive constriction,El the SGS survives at zero 
temperature even for perfect NS interfaces. In this case, 
the SGS occurs due to coherent impurity scattering of 
quasiparticles inside the proximity region, with the am- 
plitude approximately proportional to the characteristic 
length £a of this region. If we neglect the proximity cor- 
rections (m± — > 0), the SGS disappears, along with the 
excess current, aud the I-V characteristic shows perfect 
Ohmic behavior.El Thus, we conclude that in the both 
cases of resistive and transparent NS interfaces, SGS ap- 
pears due to the normal electron backscattering in the 
proximity region. This formally corresponds to the fi- 
nite value of the renormalized Andreev resistance of the 
interface. 



VIII. ROLE OF INELASTIC SCATTERING 

Since the relative SGS amplitude increases along with 
the NS barrier strength (though the current itself de- 
creases), one might expect systems with high-resistive 
interfaces to be more favorable for the experimental ob- 
servation of SGS. However, there exists an intrinsic re- 
striction for the effect: to provide strong noncquilibrium 
of the subgap quasiparticles inherent to MAR, the time 
Td of their diffusion through the whole MAR staircase, 
from —A to A, must be smaller than the inelastic re- 
laxation time r e . The value of Td can be estimated as 
the time for diffusion over an effective length (2A/eV)d. 
At low interface resistance, Rsn <C Rn, the diffusion 
rate is basically determined by the impurity scattering: 
Td{eV) ~ (2A/eV) 2 d 2 /V. For high interface resistance, 



Rsn S> Rn, the large Andreev resistance Ra ~ Rsn/W 
present bottleneck which renormalizcs the diffusion coef- 
ficient in Td by a small factor WRn/Rsn- In this way, 
the level of noncquilibrium of the subgap quasiparticles 
is controlled by the parameter 



W £ 



R 



N 



r d {2A) R N + Rsn/W 



(66) 



which must be large enough to allow observation of at 
least a few subharmonics in I(V), i.e., the condition 
W s > 1 determines the lower boundary for the barrier 
transparency. An analogous estimate for the inelastic 
parameter, with the barrier resistance Rt S> Rn substi- 
tuting for the Andreev resistance Rsn/W in Eq. J6^), 
was obtained in Ref. [?] for an SNINS structure with per- 
fect NS interfaces. In this case, the tunnel barrier does 
not affect the Andreev reflection but produces renormal- 
ization of the diffusion coefficient, T> — > T>Rn /Rt- 

At eU/2A < WT 1/2 , when r £ < T d (eV), the normal 
channel may be considered as a reservoir with a certain 
effective temperature (depending on the details of inelas- 
tic scattering which are beyond our model approach) , and 
the I-V curve becomes structureless. At small W e , the 
SNS junction behaves at two SN junctions connected in 
series through the equilibrium normal reservoir. 



IX. SUMMARY 

We have developed a consistent theory of incoher- 
ent MAR in long diffusive SNS junctions with arbitrary 
transparency of the SN interfaces. We formulated a 
circuit representation for the incoherent MAR, which 
may bje-,considered as an extension of Nazarov's circuit 
theorytil to a system of voltage biased superconducting 
terminals connected by normal wires in the absence of su- 
percurrent. We constructed an equivalent MAR network 
which includes a new resistive element, "Andreev resis- 
tor" , which provides electron-hole conversion at the SN 
interfaces. Separate Kirchhoff rules are formulated for 
electron (hole) population numbers and diffusive flows. 
Within this approach, the electron and hole population 
numbers are considered as potential nodes connected 
through the tunnel and Andreev resistors with a dis- 
tributed voltage source - the equilibrium Fermi reservoirs 
in the superconducting terminals. 

The theory was applied to calculation of the I-V char- 
acteristics. The subgap current decreases step- wise with 
decreasing applied voltage in junctions with resistive in- 
terfaces, while in junctions with transparent interfaces 
an appreciable SGS appears only on the differential resis- 
tance. In all cases, dV/dl exhibits sharp structures whose 
maximum slopes correspond to the gap subharmonics, 
eV n — 2A/n. The amplitude and the shape of SGS os- 
cillations strongly depend on the interface/normal-metal 
resistance ratio r = Rsn/Rn and reveal a noticeable 
parity effect: difference of the amplitudes of the even and 
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odd structures. This effect is specific for diffusive junc- 
tions: it comes from the strongly enhanced probability 
of Andreev reflection at zero energy. Inelastic scattering 
results in smearing of the SGS, which disappears at small 
applied voltage. 

Our theory of incoherent MAR is valid as soon as the 
applied voltage eV is much larger than the Thouless en- 
ergy: in this case, one may neglect the overlap of the 
proximity regions near the NS interfaces. In the op- 
posite case, eV < i?Th <C A, the problem of the low- 
energy part of the effective circuit should be considered 
more carefully, by taking into account the interference 
between the proximity regions. Aspects of this ac Joseph- 
son regime have been considered in Ref. [2^ in terms of 
adiabatic oscillations of the quasiparticlc spectrum and 
distribution fpactions, which produce nonequilibrium ac 
supercurrent.O Within our approach, this effect will in- 
troduce an effective boundary condition for the probabil- 
ity currents at small energy which must be included into 
the circuit representation of MAR in energy space. 

It is useful to discuss the effect of dephasing on MAR 
which has not been sufficiently investigated yet. In the 
present case of diffusive junction, the dephasing is mod- 
eled by an imaginary addition iT to the energy, E — > 
E + iT. It is easy to see that this model leads to more 
pronounced SGS. Indeed, as it follows from Eq. (Q), the 
effect of large dephasing rate T is similar to the effect of 
an opaque interface barrier (the transmissivity parameter 
W becomes effectively small). This is easily understood 
since the dephasing suppresses the anomalous Green's 
function within the normal metal similar to the effect of 
the interface barrier. Thus, even for highly transmissive 
interfaces, the I-V curves for large T ^> E^h become 
similar to the ones in the case of low-transparent inter- 
faces (see Fig. ^|), with deficit current and pronounced 
SGS. This conclusion is also supported by direct numer- 
ical calculation. It is of interest that the parity effect in 
SGS almost disappears due to the cutoff of the long-range 
proximity effect at small energies. In such case, the inco- 
herent MAR regime persists at arbitrarily small voltages 
and our theory is valid until the quasiparticle relaxation 
will affect MAR regime as described in Sec. VIII. 



RsnGi 



EQ(E-A) WA 2 Q[±(E - A) 




E 2 - A 2 

WA ( 
1 



m± 



CaW 2 A 2 



VA 2 - E 2 



2e) ^ 



(A2a) 



R N 2d\E 2 - A 2 

+Q(E - A) - 0(A - E)] 



(A2b) 



where 0(x) is the Heaviside function and E is assumed 
for brevity to be positive. From Eq. (A2a), we obtain 
approximations for the tunnel and Andreev resistances: 



R T (E)^R SN (1-A 2 /E 2 ) 



(A3a) 



Ra{E) 



2Rsn\E 



WA 



1 + WQ(A - E) \2EA/{E 2 - A 2 )| 1/21 




(A3b) 



In the vicinity of the gap edges, \E — A\ < W 2 , and 
at small energies, E < W 2 , where 9n{E) in Eq. (jAJ) 
diverges, the following approxi mate solutions of Eq. (|7) 
have to be used instead of Eq. (Al): 



exp(6W) = u 2 (t), t = 2\E- A\/AW 2 



(A4) 



at \E — A\ <C A, where u(t) is the solution of a cubic 
equation u 3 - 

sinh -y = - exp (- Arcsinh v / -iE/2W 2 A^j (A5) 

at \E\ <C A. The asymptotics of m± = ±i?jv(^A/^)A i ± 
and G± near these "dangerous" points, are given by 



RsnG+ = (Vst-^ 6 /2W^ Q(E — A), 



(A6a) 
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APPENDIX: 



The analytical expressions for bare conductivities G± 
and proximity corrections m± can be obtained in the case 
of low-transparent NS interface, W <ti 1, by making use 
of a perturbative solution of Eq. (|^) : 



6 N {E) = W y/iA/E sinh 9s{E), 



(Al) 



R SN G- = (t- 5 / 6 /4W)\V3e{E-A) + e(A-E)] , (A6b) 



// + =V2/3tV3 ^2+V3+V2-V3j , n-=\/2 (A6c) 

at t < 1. At E < W 2 , 

RsnG- = 1, M _ = (V2-l)^/A/E, (A7) 

For perfect SN interface, G^ 1 = 0, the asymptotics of 
I f -(E ) at E - A < A and E < A are given by Eqs. 
( |A6c ) , ( |A7D respectively, whereas m + (E) diverges at the 
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gap edge as [A/2(E — A)] 1 / 4 . Several examples of these 
dependencies calculated numerically are presented in Fig. 
(see discussion in Sec. IV). 



1 T. M. Klapwijk, G. E. Blonder, and M. Tinkham, Physica 
B+C 109-110, 1657 (1982). 

2 M. Octavio, M. Tinkham, G. E. Blonder, and T. M. Klap- 
wijk, Phys. Rev. B 27, 6739 (1983); K. Flensberg, J. Bin- 
dslev Hansen, and M. Octavio, ibid. 38, 8707 (1988). 

3 G. B. Arnold, J. Low Temp. Phys. 68, 1 (1987); U. Gun- 
senheimer, and A. D. Zaikin, Phys. Rev. B 50, 6317 (1994). 

4 E. N. Bratus', V. S. Shumeiko, and G. Wendin, Phys. Rev. 
Lett. 74, 2110 (1995); D. Averin, and A. Bardas, ibid. 75, 
1831 (1995); J. C. Cuevas, A. Martin-Rodero, and A. Levy 
Yeyati, Phys. Rev. B 54, 7366 (1996). 

5 A. Bardas, and D. V. Averin, Phys. Rev. B 56, 8518 (1997); 
A. V. Zaitsev, and D. V. Averin, Phys. Rev. Lett. 80, 3602 
(1998). 

6 A. V. Zaitsev, JETP Lett. 51, 41 (1990); Physica C 185- 
189, 2539 (1991). 

7 E. V. Bezuglyi, E. N. Bratus', V. S. Shumeiko, and G. 
Wendin, Phys. Rev. Lett. 83, 2050 (1999). 

8 A. F. Volkov, and T. M. Klapwijk, Phys. Lett. A 168, 217 

(1992) . 

9 W. M. van Huffelen, T. M. Klapwijk, D. R. Heslinga, M. 
J. de Boer, and N. van der Post, Phys. Rev. B 47, 5170 

(1993) ; M. Kuhlmann, U. Zimmerman, D. Dikin, S. Abens, 
K. Keck, and V. M. Dmitriev, Z. Phys. B 96, 13 (1994); 
U. Zimmerman, S. Abens, D. Dikin, K. Keck, and V. M. 
Dmitriev, ibid. 97, 59 (1995); R. Taboryski, T. Clausen, J. 
Bindslev Hansen, J. L. Skov, J. Kutchinsky, C. B. S0rensen, 
and P. E. Lindelof, Appl. Phys. Lett. 69, 656 (1996); J. 
Kutchinsky, J. Taboryski, C. B. S0rensen, A. Kristensen, 
P. E. Lindelof, J. Binslev Hansen, C. Schelde Jacobsen, and 
J. L. Skov, Phys. Rev. Lett., 78, 931 (1997); J. Kutchin- 
sky, R. Taboryski, O. Kuhn, C. B. S0rensen, P. E. Lindelof, 

A. Kristensen, J. Bindslev Hansen, C. Schelde Jacobsen, 
and J. L. Skov, Phys. Rev. B 56, 2932 (1997); R. Tabo- 
ryski, J. Kutchinsky, J. Bindslev Hansen, M. Wildt, C. 

B. S0rensen, and P. E. Lindelof, Superlattices Microstruct. 
25, 829 (1999); A. Frydman, and R. C. Dynes, Phys. Rev. 
B 59, 8432 (1999); Z. D. Kvon, T. I. Baturina, R. A. Do- 
naton, M. R. Baklanov, K. M aex, E. B. Olshane tsky, A. 
E. Plotnikov, and J. C. Portal, |cond-mat/9907247| (unpub- 
lished) . 

10 G. Johansson, G. Wendin, K. N. Bratus', and V. S. Shu- 
meiko, Superlattices Microstruct. 25, 905 (1999). 

11 Yu. V. Nazarov, Phys. Rev. Lett. 73, 1420 (1994); Super- 
lattices Microstruct. 25, 1221 (1999). 

12 A. I. Larkin, and Yu. N. Ovchinnikov, Sov. Phys. JETP 
41, 960 (1975); 46, 155 (1977). 

13 M. Yu. Kupriyanov, and V. F. Lukichev, Sov. Phys. JETP 
67, 89 (1988). 

14 G. E. Blonder, M. Tinkham, and T. M. Klapwijk, Phys. 
Rev. B 25, 4515 (1982). 

15 C. J. Lambert, R. Raimondi, V. Sweeney, and A. F. Volkov, 
Phys. Rev. B 55, 6015 (1997). 



16 A. D. Zaikin, and G. F. Zharkov, Sov. J. Low Temp. Phys 
7, 375 (1981). 

17 E. V. Bezuglyi, E. N. Bratus', and V. P. Galaiko, Low 
Temp. Phys. 25, 167 (1999). 

18 S. N. Artemenko, A. F. Volkov, and A. V. Zaitsev, Sov. 
Phys. JETP 49, 924 (1979). 

19 Here we apply the term "Andreev resistance" (and further 
"Andreev conductance" ) to denote particular values of Ra 
(or Ga). This must be distinguished from the term "An- 
dreev conductance" widely used in literature for physical 
conductance of the junction at the subgap applied voltages. 

20 A. F. Volkov, A. V. Zaitsev, and T. M. Klapwijk, Physica 
C 210, 21 (1993). 

21 In the terms of the circuit theory, this means that the "volt- 
age drop" between the electron and hole nodes is directed 
against the probability current flowing through the An- 
dreev resistor. 

22 Yu. V. Nazarov, and T. H. Stoof, Phys. Rev. Lett. 76, 823 
(1996); T. H. Stoof, and Yu. V. Nazarov, Phys. Rev. B 53, 
14496 (1996); A. F. Volkov, N. Allsopp, and C. J. Lambert, 
J. Phys.: Cond. Matter 8, L45 (1996). 

23 S. V. Lempitskii, Sov. Phys. JETP 58, 624 (1983); F. Zhou, 
and B. Spivak, JETP Lett., 65, 369 (1997); N. Argaman, 
Superlattices Microstruct. 25, 862 (1999). 

24 K. W. Lehnert, N. Argaman, H.-R. Blank, K. C. Wong, S. 
J. Allen, E. L. Hu, and H. Kroemer, Phys. Rev. Lett. 82, 
1265 (1999). 



14 



